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Abstract 

Gravitational waves (GWs) from the inspiral of compact remnants (CRs) into massive black holes (MBHs) will 
be observable to cosmological distances. While a CR spirals in, 2-body scattering by field stars may cause it 
to fall into the MBH before reaching a short period orbit that would give an observable signal. As a result, 
only CRs very near (~ 0.01 pc) the MBH can spiral in successfully. In a multi-mass stellar population, the 
heaviest objects sink to the center, where they are more likely to slowly spiral into the MBH without being 
swallowed prematurely. We study how mass-segregation modifies the stellar distribution and the rate of GW 
events. We find that the inspiral rate per galaxy for white dwarfs is 30 Gyr -1 , for neutron stars 6 Gyr , and for 
stellar black holes (SBHs) 250 Gyr -1 . The high rate for SBHs is due to their extremely steep density profile, 
jibhW oc r~ 2 . The GW detection rate will be dominated by SBHs. 

Subject headings: black hole physics — stellar dynamics — gravitational waves — Galaxy: center 



1. INTRODUCTION 

Massive Black Holes (MBHs) with masses M, < 5 x 
1O 6 M have Schwarzschild radii r s = 2GM,/c 2 , such that 
a test mass orbiting at a few rs emits gravitational waves 
(GWs) with frequencies 10 -4 Hz < v < 1 Hz, detectable 
by the planned space based Laser Interferometer Space An- 
tenna 2 (LISA). Main sequence (MS) stars with mass M* 
and radius will be disrupted at the tidal radius r t = 
(M./M*) 1 ' 3 i2* > rs and are therefore unlikely to be sources 
of observabl e GWs (our own Galactic center may be an excep- 
tion, Freitag 2003). Compact remnants (CRs) such as white 
dwarfs (WDs), neutron stars (NSs) and stellar black holes 
(SBHs) have tidal radii r t < rs and can emit GWs that are 
observable to cosmological distances. The inspiral of a CR 
into a MBH ("extreme mass ratio inspiral sources" [EMRIs]) 
is among the main targets of LISA . 

The event rate of EM RIs ha s been e stimated by nu mer- 
ous a uthors llHils & Bender 1995; Sigurdsson & Ree; 
1 19971 iMiralda-Escude & Gould 1200(1 IFreitad 1 2001 1 
Ivanovl l2002t IFreitad 120031: lAlexander & Hopmanl EPOS 
Hopman & Alexander 2005, 2006) but remains rather uncer- 
tain, in part because of the slow nature of the inspiral process, 
which occurs on many dynamical times. This makes the 
inspiraling star very susceptible to scattering by other stars, 
which can change the orbital parameters. The formalism 
for inspiral rates is similar to that for prompt consumption 
of stars teahcall& .Wollll977t iLightman & Shamrolll977t 
[ Frank & Reesl 119761 ICohn & Kulsrud 119781 ISver & UlmeJ 
1999; Magorrian & Tremaine 1999), but there are some 
important differences because the process is much slower. 

The picture can be understood as follows: Let t r be the 
relaxation time of a star with negative energy E (hereafter 
"energy"; E > for bound stars) and specific angular mo- 
mentum J (hereafter "angular momentum"). The relaxation 
time is the time-scale for a change of energy of order E, or 
a change in angular momentum of order J C (E), the circular 
angular momentum. The change in J of a star per orbital pe- 
riod P is A J = (P /tr) 1 / 2 J c . The time-scale for a change of 

1 The William Z. & Eda Bess Novick career development chair 
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order J is t j = (J/ J c ) 2 t r . In particular, the time-scale for a 
change in J by the order of the loss-cone, determined by the 
angular momentum of the last stable orbit Jlso = 4GM m /c, 
is ti c = {Jlso/ Jc) 2 t r - Inspiral due to dissipation by GW 
emission happens on a time-scale to(E, J), which for highly 
eccentric orbits has a very strong angular momentum depen- 
dence, fo(J) oc J 7 . If ti c "C to(E, J — > Jlso)> the angular 
momentum will be modified even if the star has J > Jlso- 
As a result it is very likely that the star will be scattered into 
the loss-cone (or away from it, to an orbit where energy dissi- 
pation is very weak). Such CRs will eventually be consumed 
by the MBH and add to its mass, but they will not be observ- 
able as GW emitters (GW bursts in our own G C may form an 
exception [Rubbo, Holley-Bockelmann & Finn 2006]). 

The approximate condition to(E,J — ► Jlso) < Uc{E) 
translates into a minimal energy or maximal semi-major axis 
oqw a CR must have in order to spiral in and become a 
LISA s ource ("successful inspiral"); Hopman & Alexander 
d2005t hereafter HA05) estimate that for a MBH of M. = 
3 x 10 6 M Q , acw ~ 0.01 pc: nearly all CRs with a > a G w 
are promptly captured or deflected without giving an observ- 
able signal, while nearly all stars with a <C acw do spiral in 
successfully. 

The fact that the distribution of CRs near MBHs is crucial 
to the observational outcome, implies that mass-segregation 
is likely to play a very important role for EMRIs. Mass- 
segregation is a manifestation of dynamical friction. It drives 
the heaviest objects to the center, so their concentration within 
acw increases, and drives the lightest stars to larger radii, 
so that they are relatively rare within oqw- The importance 
of mass-segregation on inspiral proc esses was dramaticall 
demo nstrated in A^-body simulations (Baumgardt et al. 200' 
1 20051) of tidal captur e of MS stars jAlexander & Morrisl2003l 
Homri an et all 12004ft . iBaumgardt et all (120051) studied tidal 
capture by a ~ 10 3 M Q black hole in a young stellar cluster 
with MS masses up to ~ 100M Q . In spite of the fact that 
massive stars are scarce, captured stars typically had masses 
- 2OM . 

In this Letter we study the implications of mass-segregation 
on the EMRI rate. 

2. MODEL 
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Our model is based on Bahc all & Woi3jl976lH977h . Here 
we briefly recapitulate the main assumptions, and discuss our 
treatment of GW capture. A detailed discussion of our model 
can be found in HA05 and Hopman & Alexander (2006). 

2.1. Dynamics 

The MBH dominates the dynamics of stars within its 
"Bondi radius", or radius of influence, = GM % /a1, where 
o\t is the velocity dispersion of a typical star of mass M* <C 
M, (assumed of Solar type), which we will use to scale our 
expressions. Orbits are assumed to be Keplerian within r/,. 
Each species with mass M is described by a distribution func- 
tion (DF) in energy space /m(E). 

We define a dimensionless time r = t/Th in terms of the 
relaxation time at the radius of influence 

3(2^/3/M*) 3 / 2 



T h = 



(1) 



327r 2 G 2 M 2 n* In A' 

where is the number density at for the typical star 
M*, (5 = Aha 2 , and A = M./M*. Introducing the dimen- 
sionless energy x = (M+/M)(E/ (3) and the dimensionless 

DF g M (x) = [(27r/3/M^) 3 / 2 n- 1 ]/M(^), t he Fok ker-Planck 
equation in energy space is (Bahcall & Wolf 1 19771 Eg. [26]) 



dg M (x,r) 



r 5/2 



8 



}m(x) - Rm(x). 



(2) 



dr dx 
We also write equation in a logarithmic form suitable for 
numerical integration (see appendix). The spatial number 
density n m (r) of stars is related to the DF by 

n M (r) = 2ir~ 1/2 n* / dxg M {x)[r h /r - x] ' . (3) 



We fit our numerical results by power-laws um{t) oc r ~ aM . 

In expression 0, Qm(x) is the (dimensionless) rate at 
which stars flow to energies larger than x, 



^ M M> 

M' 



9m{x) 



dx' [max (a;, x')] 3 ^ 2 x 
dg M '{x') 



dx' 



M' n dg M {x) 
-M9M>(x)-^- 



(4) 



The dimensional stellar current is related to Qm by 

I M (E,t) = I Q M (x,t), where 



h = 



8tt 2 , (GM*) 2 liiAr!.* 



(5) 



3_\/2_ (J.. 

(Bahc all & Wolll976tlHopman & Alexandeil2006l) . 

The last term in equation (j2ji represents losses of stars due to 
loss-cone effects (both prompt infall and inspiral) in J-space. 
The sink term in the diffusive regime for the loss-cone is 



Rm(x) 



9m {x) 



(x) ln[J c (x)/ Jlso] ' 



(6) 



where J c (x)/ J LSO = (l/4V2)(c/er*)£~ 1/2 - The full-loss 
cone regime, x < 10, d oes not contribute to the GW event rate 
(Alex ander & Hopman 120031: HA05; Hopman & Alexander 
2006). In our calculations we neglect the sink term in the full 
loss-cone regime by setting Rm — > for x < 10. In Eq. (|6) 
the dimensionless local relaxation time r r (x) is 



T r (x) 



Ml 



Em9m(x)M* 



(7) 



independent of the stellar mass. 

Let the number of stars accreted to the MBH before giv- 
ing an observable GW signal be N p (x), and the number 
of those that spiral in successfully and do give a signal 
Ni(x). The steady state result for r r (x) is used to deter- 
mine the probability for inspiral Sm(x) — Ni(x)/[Ni(x) + 
N p (x)] by Monte Carlo simulations (HA05) as follows. At 
every orbit, a star of initial energy E and initially large 
J makes a step in J of order A J = [P(x)/<: T .(x)] 1 / 2 J c 
with random sign because of scattering, and loses energy 
AEgw = ( 857r/24576)(M/M.)Mc 2 (J/J T ,so)~ 7 to GWs 
(Petersl [T964) . This is repeated many times and the outcome 
is recorded. The total rate of successful inspirals for species 
M is then given by T a/ = h J°° dxSM(x)x~ b / 2 Rm{x). It 
is convenient to express the capture rate in terms of the semi- 
major axis a = r/,./2x of the stars, 



r Af (< a) 



2V2I0 

3/2 



daa 1/2 S M (a)R M (a) . (8) 



2.2. Boundary conditions and model parameters 

Equation (0 has inner and outer boundary conditions. At 
some large energy ijj the DF vanishes, g(x > Xd) = 0. Since 
the EMRI rate is dominated by the largest distance where suc- 
cessful capture is possible, the exactly value of xd is not im- 
portant. Here we used xd = 10 4 , which is approximately the 
energy-scale where the inspiral time becomes smaller than a 
Hubble time even for a circular orbit. We assume that the 
MBH mass is M. = 3 x 10 6 M Q , representative of a typical 
LISA source. Our value of xd would approximately corre- 
spond to a distance scale of ~ 10 -4 pc from the MBH. 

The second boundary c ondition is given at x = 0: follow- 
ing Bahcai i"& Wolll ( 119771) . we assume that for x < the 
stars have a Maxwellian velocity DF with equal temperature 
(Pm = Mctjtf — (3), and with different population number 
fractions for different species Cm, 



g M (x) = C M exp(Mx/M*); (x < 0). 



(9) 



We consider four populations of stars. One species consists of 
main sequence stars, assumed here to be of Solar mass. MSs 
do not contribute to the GW inspiral rate since they are tidally 
disrupted before spiraling in, but they do contribute dynam- 
ically and they dominate both in number and in total mass 
at the radius of influence. The other three populations con- 
sists of WDs (M WD = 0.6M©), NSs (M NS = 1.4M ) and 
SBHs (A/bh = IOMq). The number fraction ratios of the 
four populations at x = are Cms : Cwd : Cns : Cbh = 
1 : 0.1 : 0. 01 : 10~ 3 , typica l for continuously star forming 
populations llAlexanderll2005i) . We also adopt for our model 
the Galactic center value s a+ = TSkm s" 1 (77, = 2pc) and 
n* = 4 x 10 4 pc" 3 (lGenzel etalJ20 03) 3 . The model param- 
eters are summarized in table Q. 

3. RESULTS 

We integrated Eq. (0 until steady state is obtained, af- 
ter time t < 1. In figure Q we show the resulting densi- 
ties for the different species. The DF of the SHBs is much 
steeper than that of the other types (q;bh = 2.0), and at 



3 The Galactic MBH obeys the M. - a r elation 1 Ferrarese & Merr itl200(l 
IGebhard t et al. 2000; Tremaine et al. 2002), so that these values may also be 
representative of other MBHs. 
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FIG. 1 . — Steady state number densities njy/M- The DF of SBHs (dashed- 
dotted) is much steeper than that of the other three stars (obh ~ 2). SBHs 
dominate NSs (dotted) by number nearly everywhere, and they even dominate 
WDs (red dashed) within ~ 0.01 pc, in spite of their low number fraction at 
Th> C'bh/C'wd = 1%. MS stars (solid) dominate by number everywhere. 
For the other three species we found {omSi "WD, cvns} = {1-4, 1.4, 1.5}. 
The relaxation time (rising solid line) grows approximately as t r cxr° . 



Table 1 

Model parameters and GW rates 



Star 


M 


Cm 


O-M 


7V A f(<0.01pc) JV M (<0.1pc) 


fl GW Vm 




[Mq] 










[mpc] [Gyr" 1 ] 


MS 


1 


1.0 


1.4 


10 3 


3 x 10 4 




WD 


0.6 


0.1 


1.4 


80 


2.7 x 10 3 


4 30 


NS 


1.4 


0.01 


1.5 


12 


374 


5 6 


SBH 


10 


10~ 3 


2.0 


150 


1.8 x 10 3 


13 250 



r ~ 0.01 pc the number density of SBHs becomes compa- 
rable to that of the WDs. MS stars dominate everywhere by 
number, althou gh we did not take into account stellar colli- 
sions (Freitag & Benz 2002, 2005) which could deplete the 
MSs close to the MBH. SBHs also determine the functional 
behavior of t r oc r p , where pwaBH - 3/2 « 0.5. 

Throughout most of the cusp ckbh ^ 2, but near ijj the DF 
flattens, as required fo r the integrals in equa tion to con- 
verge at high energies (Bahcall & Wolf 1 1 9771) . Large slopes 
at intermediate energies are allowed by these equations, and 
arise when a population of massive objects with a low num- 
ber density exists, as is the case in our model. At low energies 
the massive objects sink effectively to the center by dynam- 
ical friction. At high energies the massive objects dominate 
the dynamics, decouple from the lighter objects, and form an 
a = 7/4 "mini-cusp". This process is reminiscent of the 
Spitzer instability in glo bular clusters, where SBHs decou- 
ple fr om the other stars (Spitzer & HartlH971l iKhalisi et alJ 
2006). 

The probability for inspiral Sm(o) is shown in figure (0. 
Since the SBHs are more massive they lose energy to GWs 
at a higher rate than the other species and can spiral in from 
larger distances. 

In figure we show the cumulative rates of successful 
inspiral (Eq. [8jl for all CRs as a function of distance from the 
MBH. We summarize some results in table Q, where we also 
give the enclosed number of stars Nm(< a) within a. 

4. SUMMARY AND DISCUSSION 
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FIG. 2. — Probability of successful inspiral for a consumed star as a function 
of distance from the MBH for the CRs. BHs (dotted-dashed) can successfully 
spiral in from further distances than WDs (solid) and NSs (dashed) due to 
their higher masses. 
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FIG. 3. — Cumulative rates as a function of a are dominated by BHs (dotted- 
dashed); WDs (solid) and NSs (dashed) contribute significantly less. 



In our model, SBHs dominate the EMRI rate, in spite of 
their small number density at rn- The combination of a very 
steep cusp (q:bh ~ 2.0) and a larger oqw due to their larger 
mass leads to Tbh > (Fwd, Tns) P er galaxy. We also note 
that the amplitude of the GWs is proportional to the stellar 
mass, so that the distance at which these objects can be ob- 
served is~ 10 times larger than that for WDs and NSs. Thus, 
SBHs will dominate the cosmic detection rate. 

It is instructive to compare the EMRI rates we obtain here to 
those obtained by HA05, where mass-segregation was not ex- 
plicitly included. For SBHs, Tbh is larger by a factor ~ 50. 
Part of the difference is that we assume here a larger total 
number of SBHs within the cusp (by a factor ~ 6; we nor- 
malized the SBH number fraction at r/j to be Cbh = 10~ 3 , 
while HA05 assumed that the enclosed fraction of SHBs is 
10~ 3 ). More importantly, the steeper cusp leads to a higher 
capture rate (by a factor ~ 9, HA05, eq. [32]; HA05 assumed 
«bh = 175). The BH c usp is much s teeper than any of the 
cases studied by Bahcall & Wolf ( 1977), and in particular it is 
steepe r than the cusp o f a single mass population, a = 7/4 
(Bah call &WoHl976l) . 

The rates Twd and Tns are also somewhat larger than those 
found by HA05. Here the difference originates mainly in the 
behavior of t r : For WDs and NSs, t r was assumed to be 
constant by HA05, as appropriate for a single mass popula- 
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tion with a = 3/2. However, the interaction between SBHs 
and the other CRs leads to a decrease in t r towards the MBH 
(Fig. Using the analytical expressions by HA05, it can 
be shown that if nA/(r) oc r~ 3 / 2 , and t r oc r p , the suc- 
cessful inspiral rate is enhanced by (d c /rf- l )~ 3p ^ 3 ~ 2p ^ ~ 10 
(for p = 0.5) relative to the t r — const, case, where d c = 
[(85/3072)V^M;(M/M.)^] 2 / 3 , see HA05 eq. (29). 

The EMRI rates w e found here are promising for 
the L ISA detection rate ( Bar ack & Cutlerl I2004b1: iGair etafl 
2004), in spite of the fact that more sources also imply a 
stronger background noise (Barac k~& Cutlerl2004al) . 

We neglected here the effect of resonant relaxation 
llRauch & Tremainell996URauch & Ingallsll9 98). which can 
increase the EMRI rate by up to an order of magnitude 
(Hopman & Alexander 2006). A multi-mass analysis of RR 
has yet to be performed. In addition to direct capture of 
CRs, EMRI ca n occur follo wing the formation of SBHs in 
accr etion disks tLevinl2 003). binary disruptions ( M iller et alJ 
120051) and tidal capture follo wed by a super nova exp losion of 
the captured star (Hopman & Portegies Zwart 2005). These 
other mechanisms lead to low eccentricity signals, whereas 



direct capture leads to high eccentricities (HA05). 

Our estimate of the number fraction of unbound SBHs is 
somewhat uncertain, in part because we neglected dynam- 
ical effects for unbound stars. We note that our estimate 
-Mbh(< Pc) ~ 1.6 x 10 4 is c onsistent with calculations by 
Miralda-Escude & Gould fcOOOl) . who found 7V"bh(< pc) ~ 
2.5 x 10 4 . Our Galactic Center contains a M BH of mass 
comparable to the MBH mas s considered here JGhez et alJ 
120051 Eisenhau er et alJ 12005). Observational effects of a 
clus ter of SBHs near the Galactic MBH includ e microlens- 
ing (IChaname et afl bOOl). X-ray emission iPessah & Meli al 
2003), capture of massiv e stars by an exchange interaction 
(Alexander & Livio 2004) and deviations from Keplerian mo- 
tion of luminous stars (Mouawadet al. 2005). Such effects 
could in principle be used to constrain the predicted densities. 
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APPENDIX 



Because of the large range of energies, the natural way to integrate the Fokker-Planck equation is to divide the energy range 
into equal logarithmic intervals. For convenience we give here the equations in terms of the logarithmic distance variable z = 
ln(l + x/X). The Fokker-Planck equation without sink terms is then written as 



dg M (z,T) _ M_ n 5/2- z 



E 



M 



M' d r 
1 ~Ml~d~z 



(e z - l)- 3 l 2 g M {z)[g M ,{z) - g M ,(-oo)] + g M {z) f ° dz'(e z ' - i)-3/a gg*dg) 



M' 
~M 



dz 



, , M' _ z dg M (z) 
dze 9 M>(z)-—e 



dz' 

dz'e z '(e z ' -ir 3 / 2 g M >(z'). 



(1) 



The logarithmic expressions for the sink term (eq. |6) can be included directly by replacing x — » X(e z — 1). 
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